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While of paramount importance in material science, the dynamics of cracks still lacks a complete 
physical explanation. The transition from their slow creep behavior to a fast propagation regime is a 
notable key, as it leads to full material failure if the size of a fast avalanche reaches that of the system. 
We here show that a simple thermodynamics approach can actually account for such complex crack 
dynamics, and in particular for the non-monotonic force-velocity curves commonly observed in 
mechanical tests on various materials. We consider a thermally activated failure process that is 
coupled with the production and the diffusion of heat at the fracture tip. In this framework, the 
rise in temperature only affects the sub-critical crack dynamics and not the mechanical properties 
of the material. We show that this description can quantitatively reproduce the rupture of two 
different polymeric materials (namely, the mode I opening of polymethylmethacrylate (PMMA) 
plates, and the peeling of pressure sensitive adhesive (PSA) tapes), from the very slow to the very 
fast fracturing regimes, over seven to nine decades of crack propagation velocities. In particular, the 
fastest regime is obtained with an increase of temperature of thousands of kelvins, on the molecular 
scale around the crack tip. Although surprising, such an extreme temperature is actually consistent 
with different experimental observations that accompany the fast propagation of cracks, namely, 
fractoluminescence (i.e., the emission of visible light during rupture) and a complex morphology of 


post-mortem fracture surfaces, which could be due to the sublimation of bubbles. 


I. INTRODUCTION 


The rupture of solids is often described by empirical 
observations rather than by fully understood physical 
models. One of the earliest formalisms is that by Grif- 
fith in 1921[1]: the propagation of cracks is described 
as a threshold phenomenon, only obtained when loading 
their encompassing matrix above a critical fracture en- 
ergy. To the first order, this view matches the behavior 
of brittle bodies, which suddenly snap passed a certain 
elastic deformation. Analytical models of cracks propa- 
gating in lattices suggested [2—5] that such an instability 
arises from the discrete nature of matter at the atomic 
scale. Indeed, these models revealed a minimum prop- 
agation velocity, comparable to that of the mechanical 
waves in the considered material, above which the ad- 
vance of a fracture tip through the network of molecular 
bonds can be self maintained by the emission of high 
frequency phonons. There, the energy binding two lat- 
tice nodes is defined as a covalence-like barrier [6]. While 
this description [2] does not allow for slow propagation, 
it is acknowledged that a crack loaded well below the 
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fast rupture threshold is still growing, but at creeping 
rates that are orders of magnitude below that of a ‘dy- 
namic’ fracture (e.g., [7, 8]). An approach to explain such 
a creep in a way that is compatible with Griffith’s for- 
malism [1] is to consider that the fracture energy is not 
an intrinsic material property, but is instead a particular 
function of the propagation velocity (e.g., [9]). One hence 
simply obtains a lower crack speed if providing a lesser 
mechanical load. Alternatively, the creep regime is well 
modelled [7, 8, 10-16] by thermally activated sub-critical 
laws such as Arrhenius-like growth rates (e.g., [17]), and 
thermodynamics has thus emerged as a framework to de- 
scribe the slow failure. In such descriptions, that are 
sometimes referred to as ‘stress corrosion’, a variation of 
fracture energy with velocity is not particularly called for, 
as the molecular agitation allows the crack to progress at 
loads below an intrinsic rupture threshold. 

In practice, and depending on the material being bro- 
ken, both the slow and the fast propagation regimes can 
be observed for a same range of applied loads[18, 19]. 
A hysteresis holds and the growth rate of a fracture is 
then depending on the actual mechanical history, rather 
than only on the instantaneous mechanical load. Maugis 
and Barquins [20, 21] early suggested that the descrip- 
tion of the slow and the fast regimes, as well as that 
of the hysteresis, could be qualitatively unified by rein- 
terpreting Griffith’s criteria |1], if one could account for 


the temperature and velocity dependent viscoplasticity 
that occurs around crack tips [22, 23]. More specifically, 
Marshall et al. [18] and then Carbone and Persson [24, 25] 
proposed that the induced heat associated to such a plas- 
ticity might locally soften the matter around a crack and 
that some thermal weakening (i.e., the abrupt transition 
from slow creep to fast failure due to a thermal process) 
arises from the related reduction of the material elastic 
moduli. 

In this work, we propose a quantitative unifying model 
of the two propagation regimes that disregards such a 
softening effect, hence stating that some variations in 
the material mechanical properties are not necessarily 
required to obtain a slow-to-fast-crack transition. We 
focus instead on how the thermal dissipation, and the 
subsequent rise in tip temperature, affect the front sub- 
critical growth, as understood by statistical physics and 
an Arrhenius-like law. In some previous works, we indeed 
studied how such sub-critical laws, at fixed room temper- 
ature, well describe creep; in fibrous and polymeric ma- 
terials (namely, paper sheets and polymethylmethacry- 
late, PMMA), they notably account for the mean kinet- 
ics of slow rupture fronts under various loading condi- 
tions [14, 15, 26]. When, in addition, taking into account 
these media structure and heterogeneities in fracture en- 
ergy, such sub-critical laws also reproduce the intermit- 
tent dynamics of failure; in particular, the size distribu- 
tion of crack jumps [13] and the front roughening prop- 
erties [16]. Here we neglect any spatial variation of the 
fracture energy, but let the crack tip temperature vary 
as a function of the front velocity and of the applied me- 
chanical load. Indeed, in a previous experimental and 
theoretical study of the tearing-induced heating in paper 
sheets [27], we were able to relate the temperature field 
around moving cracks to a certain percentage of the me- 
chanical energy which gets converted into heat as the tip 
advances. More recently, this rise in temperature was fed 
back into a sub-critical growth law and showed [28] that 
one can thus obtain a dynamics model holding numer- 
ous qualitative similarities with the observed behavior 
of cracks, namely, two stable phases of propagation and 
a critical point that is similar to a brittle-ductile tran- 
sition (e.g.,[29]). Here, this model is first reintroduced 
(section IT) and then shown to quantitatively capture the 
fracturing dynamics of two different polymeric materi- 
als, over the full range of velocities (section III), namely, 
acrylic glass (PMMA) and pressure sensitive adhesives 
(PSA). In both these media, some extensive experimental 
work has been carried out by different groups to quan- 
tify the two rupture regimes (e.g., see Refs. [9, 30-35] for 
PMMA and Refs. [19, 36-38] for PSA) and our proposed 
model accounts for the experimental curves of applied 
load versus crack velocity, from the slowest (microme- 
ters per second) cracks to the fastest (hundreds of me- 
ters per second) ones. Such a match suggests that the 
growth of cracks could be sub-critical (i.e., as stated by 
the model) over a far wider velocity range than what 
is commonly accepted, that is, even at propagation ve- 


locities approaching that of mechanical waves. Indeed, 
we infer that the load threshold at which cracks typi- 
cally shift to the fast phase is actually smaller than the 
intrinsic rupture energy, as a result from the boosted 
thermal activation around the front. In particular, we 
predict that crack tips can reach thousands of degrees 
on the molecular scale (i.e., over a few atoms around 
the front), when they quickly avalanche. Although such 
high temperatures are today rarely considered, they have 
long been proposed (e.g.,Rice and Levy [39]), and we here 
discuss (section IV) how they are inline with several ob- 
servables that sometimes accompany the fast propaga- 
tion of cracks, namely, the emission of visible light at 
their tips (i.e., fractoluminescence [40—42]) and the exis- 
tence of bubbles on their postmortem surfaces, that can 
nucleate secondary rupture fronts [43, 44]. 


II. FROM THERMAL CREEPING TO 
THERMAL WEAKENING 


A. The kinetics of sub-critical rupture 


We here consider a model in which the velocity V of 
cracks is ruled by the competition, at their tips, between 
breaking and healing processes [45] (or see Ref. [46], chpt 
5.5.1). As many authors before us (e.g.,[7, 8, 46]), we 
propose that these processes are, at least in part, sub- 
critical, and are governed by some Arrhenius-type laws 
(e.g., [17], chpt. 1.8.1). The activation energies of these 
laws are thus exceeded by the thermal bath according 
to a probabilistic Boltzmann distribution [17]. The rup- 
ture activation energy can then be written as (U, — U): 
the difference between the mechanical energy U that is 
stored in the tip bond and a critical rupture energy Ue, 
at which this bond fails. The latter should typically be 
comparable to a few electronvolts, which is a standard 
value for atomic covalence (e.g., see appx. E in Ref. [6]). 
Of course, depending on the studied material, Ue could 
also be dominated by the typically weaker binding ener- 
gies of hydrogen or Van der Waals bonds, and its actual 
value may thus lie within a few orders of magnitude. In 
any case, as we are here introducing a mesoscopic law 
for the rupture dynamics (i.e., an Arrhenius growth), Ue 
should be understood as a mean material property, repre- 
sentative of the various strengths of the links that break 
along a crack course. Such a statistical definition will 
also apply to most of the parameters that we will hence- 
forward consider. Similarly to the rupture barrier, the 
activation energy to heal the atomic connections can be 
written as (Un + U). There, U, is an intrinsic repulsive 
energy barrier that two atoms need overcome to bond, in 
addition to which the thermal bath at the healing link 
also needs to compensate for the applied stretch U of the 
tip. With these considerations, the propagation velocity 
of a crack is then modelled by 


U. — U Un +U 
V = vdo exp (- EBT ) = rdoexp (- a | (1) 


A Thermal 
agitation 
—— 


of NO 


Crack tip 0 


Free 


: d, d 1 Atoms 
separation 
d A d (Breaking link) 
2 | do 
Y 
A Thermal 
agitation 
soe ae NO 


é Breaking 
link © link 


> 


2 Atoms 


separation 
(Healing link) 


FIG. 1. (Left): simplified atomic view of the breaking/healing 
site at the crack tip. (Top right): Generic tip stress ø nor- 
malised by the stress shielding factor N versus atom separa- 
tion for the active breaking link. (Bottom right): Generic tip 
stress normalised by the stress shielding factor versus atom 
separation for the active healing link. The grey areas are the 
energy release rate G. At this load, dı and d2 are the mean 
extensions of, respectively, the breaking and the healing link, 
while do is the unstressed atom separation. On the breaking 
link graph: the area below the curve for d > do is the intrinsic 
surface fracture energy Ge. The thermal agitation may over- 
come the remaining Ge — G barrier. Although the healing 
link is initially broken, an energy input is required to move 
the two particles closer to each other, due to the neighbour- 
ing unbroken links stretched at a load G. In addition, when 
the atoms separation gets smaller, the thermal agitation also 
needs to overcome a repulsive energy barrier G% (the area be- 
low the atoms separation axis in this figure) before reforming 
the bond. 


where the first term is the forward rupture velocity of the 
crack and the second one is the backward healing velocity. 
In this equation, we denote do ~ 2A the inter-atomic 
distance, kg is Boltzmann’s constant ~ 1.38 x 10773 
m? kgs? Kt, T is the absolute temperature at the crack 
tip and v is the collision frequency in the molecular bath 
(e.g., [17], chpt. 4.1). Each exponential in Eq. (1) is a 
probability term (i.e., the probability, challenged every 
1/v second, that the thermal bath exceeds one of the 
activation energies and that the crack hence advances or 
retreats by a step do). As such, these terms cannot be 
greater than 1 and, while the healing one always meets 
this condition, U > Ue corresponds to an over-critical 
propagation regime where 
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The product vdo is a maximal velocity, that we will fur- 
ther denote Vo, at which a fracture front can advance, 
when its tip atomic bonds snap each time they are chal- 
lenged and never heal. In theory [17], the frequency v 


is temperature dependent, with Vo ~ \/kgT/m where 
m is the mass of an atom or a molecule, but this de- 
pendence is small compared to that of the neighbouring 
exponential terms, so that we here neglect it. In our 
context of rupture kinetics, and more practically, it was 
notably proposed [5, 47] that such a nominal velocity Vo 
is in the order of that of the medium Rayleigh waves, as 
quicker fractures then propagate in a specific supersonic 
regime [48, 49], which is not here considered. 

In our description, U is the physical quantity that de- 
scribes the load of a crack on the microscopic level, and 
that governs most of its dynamics. However, at the lab 
scale, U is not a measurable quantity. The energetic level 
at which a crack progresses is rather characterized by the 
macroscopic energy release rate G, which is the amount 
of energy that a fracture dissipates to grow by a given 
unit of measurable area [1]. This energy dissipation may 
be of diverse nature, and is to cause a relative reduction 
in potential energy near the tip. We will denote N > 1 
the factor for this reduction, so that U ~ diG/N. More 
commonly, mechanical shielding is described with the in- 
troduction of a plastic process zone of radius € around 
the crack front, where the dissipation occurs. To follow 
this canonical framework, we define a radius € that is rel- 
ative to the length of an atom link, such that 26/do = N. 
The intensity of the mechanical shielding (i.e., the rela- 
tion between the potential energy U stored in the rup- 
turing bond and the macroscopic energy dissipation G) 
then writes as 


dG 
2€ ? 
By additionally introducing Ge = 2£U./dÿ and Gr = 
2£U),/d3, the respective equivalents in the energy release 


rate framework of U, and Up, one can re-write Eqs. (1) 
and (2) as functions of G: 
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(4) 
We have here also written T as To + AT, where Tọ is 
the absolute room temperature (~ 296 K) and AT is any 
deviation from this background value, as we will proceed 
to propose that the tip temperature can vary. 
Note finally that one could also write this relation as a 
function of the mechanical stress ø that is applied at the 
crack tip, using G(d) = i No(d')dd’, where do is the 
nominal separation of atoms in an unloaded matrix (i.e., 
at G = 0) and d is the actual atom separation at the 
crack tip. Figure 1 illustrates such a link between G and 
o and summarizes, in a simplified atomistic view, how 
the thermal bath allows to overcome the surface energy 
barriers for breaking and healing atomic bonds, Ge — G 
and Gr + G, as per Eq. (4). 


B. Heat dissipation and tip temperature rise 


In the model we have introduced, one needs to further 
account for the energy which is dissipated around the 
running tip (G), as, even if it is mechanically lost, we will 
here show that it can maintain a strong effect on the crack 
dynamics. While the energy dissipation can be of several 
forms, ranging from the emission of mechanical waves [50] 
damped in the far field, to the nucleation of defaults in 
the matrix [51] (i.e., crazing [9, 52]), we here focus on the 
release of heat around the fracture tip [27, 40]. We thus 
call ¢ the percentage of G that is converted into some lo- 
cal rise in internal energy, and hence in temperature, and 
denote | the typical size over which this process occurs. 
As the heat, released on a production zone of area ml? 
close to the tip, is to diffuse in the whole bulk, the result- 
ing temperature elevation AT can be modelled (e.g., [27]) 
by the standard diffusion equation: 


AAT) À, 
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where À is the medium’s thermal conductivity, and C is 
the volumetric heat capacity. The last term of this equa- 
tion is a source term only valid in the heat production 
zone. The support function f of this zone is 1 inside of it 
and 0 otherwise, and the thermal source term is propor- 
tional to GV, that is the dissipated power per unit of 
crack length deposited in the advancing zone. Although 
governed by Eq. (5), AT at the rupture front can approx- 
imate to far simpler expressions. It was indeed shown [27] 
that, at low propagation velocities, the temperature ele- 
vation at the centre of the heat production zone (i.e., the 
crack tip) is only governed by the diffusion skin depth 
6 = \/AT/(xC) upon the passage of the production zone 
of extension l within the time r = 1/V. For fast cracks 
however, when 6 becomes smaller than l, the generated 
heat can barely diffuse out of its source zone and AT is 
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FIG. 2. Steady thermal elevation at a crack tip for various 
propagation velocities, due to the diffusion equation (5) (plain 
plot). The approximations AT fast and AT'siow, from Eqs. (6) 
and (7), are shown for comparison (dotted plots). The axes 
are not annotated for the sake of generality. 


then constrained by l. We thus have 
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Figure 2 shows the general evolution of AT at the tip 
with V, according to Eq. (5) solved by numerically inte- 
grating the heat diffusion kernel [53]. Note that AT in 
Eq. (5) is a temperature field as shown for instance in 
the inset of Fig. 3, but we are here mainly interested in 
its value at the centre of the heat production zone (i.e., 
where the rupture process occurs). Figure 2 also shows 
how the two expressions of Eqs. (6) and (7) approximate 
for the tip temperature. 


C. Model phase behavior 


We have now derived the two constitutive equations of 
our fracture dynamics model: Eq. (4), that gives the ve- 
locity of a crack as a function of its tip temperature, and 
Eq. (5), that governs the thermal state around a progress- 
ing front. In a previous work [28], we have simultaneously 
solved these two equations and, focusing on their steady 
state, showed that they predict two stable phases for the 
propagation of cracks. These two behaviors are shown 
by the plain curve in Fig. 3, and are there labelled ‘Slow 
stable phase’ and ‘Fast stable phase’. The first one, as its 
name suggests, is a slow one, where AT stays small com- 
pared to To, such that the growth rate is mainly governed 
by the medium fracture energy G, (i.e., as indicated by 
Eq. (4)). This slow branch ceases to exist beyond a par- 
ticular load G = G,. The second phase is reached when 
the generated heat (and hence AT) significantly over- 
comes the background temperature. From the Arrhenius 
law (4), the growth rate then significantly increases, so 
that the crack is said to be thermally weakened. Note, 
in Fig. 3, how both phases coexist for a certain range of 
energy release rates: a hysteresis situation holds (e.g., be- 
tween G = 300 J m~? and G = Ga in Fig. 3). When this is 
the case, the model also predicts [28] a third phase, that 
is, by contrast, unstable and hence shall be difficult to 
be recorded experimentally. 


II. COMPARISON TO EXPERIMENTAL 
RESULTS 


Interestingly, this phase description in our model 
matches key observations of fracturing experiments. The 
abrupt transition, passed a load threshold, from slow 
cracks to fast cracks, can indeed be interpreted as a phase 
transition [28], and the usual stick-slip of fronts is a good 
indicator that some hysteresis holds in the physical laws 
that rule the rupture dynamics [21, 28]. We then proceed 
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FIG. 3. Crack velocity V as a function of the energy release rate G as predicted by Eqs. (4) and (5) (plain curve) and fitted to 
the PMMA experimental data [32]. The arrows indicate to which model parameters each part of the curve is mainly sensitive, 
and the main color scale specifies at which temperature the crack tip is modelled to be. The load G, is an avalanche threshold 
beyond which a front can only propagate quickly and Ge is the modelled microscopic energy barrier for rupture. Below the 
asymptote at (Ge — Gn)/2, fronts cannot propagate forward due to some dominating healing processes. The inset shows, for 
a given point of the curve: V ~ 100ms°! and G ~ 600Jm”, the associated modelled temperature field around the front. For 
readability, the color map is there different from the main one, and the circle corresponds to the tip of radius l, where the extra 
heat is emitted. AT of the main model curve is the value at the centre of the circle. At loads beyond G = Ge, micro-cracks 
begin to nucleate [32], as shown further in Fig 4, which shows a zoom of the fast branch. 


to test our model against two sets of experimental data, 
where both the energy release rate G and the velocity V 
of the slow creep and the fast propagation stages were 
well quantified, as we detail in the next sections. 


A. The rupture of PMMA 


First, we look into a data set acquired when breaking 
polymethylmethacrylate plates (PMMA) at room tem- 
perature (Ty = 296K). A wedge is driven into Perspex® 
bodies, resulting in cracks for which two stable (G, V) 
branches are indeed recorded [32]. These results are 
shown in Fig.3. There, the fast branch, with propaga- 
tion velocities above 100 m s`}, was reported by Scheibert 
et al. [32], and the slow creeping branch is here published 
for the first time for this given PMMA (see appendix A 
for details on how it is obtained). When forcing the rup- 
ture velocity between these two regimes (i.e., above a 
specific creep velocity of 4 cms and below ~ 100 ms‘), 
some stick-slip is observed in the dynamics of the fronts, 
as reported by Hattali et al. [54]. 

Figure 3 then compares both experimental branches with 
our proposed model. We thus pursue by detailing how 
each parameter was fitted (i.e., how the model was cali- 
brated to the data), based on asymptotic read-offs. We 
classically start by wondering how well the slow propaga- 
tion phase is represented by an Arrhenius law of constant 


temperature. In the model, this corresponds to a linear 
In(V) to G relationship that holds at low velocity, where 
In is the natural logarithm. There, AT is negligible com- 
pared to the background To and G is high enough for the 
healing terms of Eq. (4) to be secondary (i.e., the terms 
involving Gp in this equation), leading to 


dè BG. 
in(v) =e | seta | + o- PSE . (8) 


In the data, this equation shall describe the portion of 
the plot lying between 107-4 and 107? mst, and the 
slope there, approximately 0.02 m? Jt, hence constrains 
d3/(2£kpTo) and so the equivalent length for the crack 
mechanical shielding € to be in the order of 50nm. Ad- 
ditionally, the intercept of Eq. (8) with the V axis (i.e., 
the second term in brackets) links Vo and Ge. We ear- 
lier stated the former to be comparable to the medium 
Rayleigh velocity [32], 880ms in this particular poly- 
mer, so that we can deduce the rupture threshold G, to 
be about 1300Jm°?. This value, together with that of 
£, gives a fracture energy Ue = Ged /(2£) comparable to 
1eV, which is satisfyingly consistent with a covalence-like 
barrier. Next, the healing threshold G, can be inferred 
from the vertical asymptote at G = 300Jm?, below 
which healing seems to prevail as cracks do not prop- 
agate forward [45]. Equation (4) predicts this asymptote 
for G = (Ge — Gh)/2, when the healing term equals the 
breaking one, such that Ga ~ 650Jm-?. Let us now fo- 


cus on the maximum G in the slow stable phase, denoted 
Ga (for ‘avalanche’) in Fig.3, around V = 4cms!. It is 
modelled by Eq. (4) once AT is high enough compared 
to To to trigger a phase transition, which, as per Eq. (6), 
mainly depends on the A/¢ ratio. By tuning this ratio, 
and appreciating the fit (see appendix B), we have de- 
duced it to be around 0.9Js!m!K°l. As the PMMA 
conductivity, À = 0.18 Js! mr! K-t, is known [55], we can 
approximate ¢ + 20%. Note that at this particular point 
(at G = Ga), the polymer suddenly breaks (e.g., [20, 56]), 
as 0V/0G — +00 and the velocity has to jump to the fast 
regime. Consequently, Ga is often seen as a macroscopic 
critical energy release rate, which in our description is 
less than the intrinsic microscopic energy barrier (i.e., 
Ga < Ge). This difference is here directly related to the 
thermal conductivity of the medium, and the avalanche 
to a fast rupture arises when the diffusion can no longer 
cope with the crack velocity, so that heat is no longer 
efficiently diffused away from the tip. The characteristic 
size l on which this heat is generated is the only parame- 
ter that remains to be determined. As, according to the 
model, the crack needs to be hot enough to explain some 
fast fronts at low mechanical load (i.e., the slower part 
of the fast branch in Fig.3, around 100 ms), we can 
estimate the limiting factor of AT fast, Cl (see Eq. (7)). 
Matching the data set in this area (see appendix B), and 
using [55] C ~ 1.5 x 10° J Km, we have deduced } to 
be in the nanometer range. This magnitude happens to 
be in the same order as the earlier derived €. We thus 
predict that most of the induced molecular agitation is 
introduced on the closest atoms around the crack tip, 
which coincides with the length scale for the energetic 
shielding of the tip. Noteworthily, such a nanometer scale 
appears to be close to the typical entanglement scale of 
polymers [57] (i.e., the density of polymeric chains cross- 
ing points in the matrix). 

To quantify how well the model accounts for the experi- 
mental data, we computed, for each data point, the rel- 
ative orthogonal distance £q to the model, that is 


2 2 
cl Gas¥e) = minm |1- Se] à [x — ln |” 
d 


where the subscript d stands for ‘data’ and m for ‘model’. 
We are thus looking at a relative fit mismatch along the 
G axis and a relative fit mismatch, in order of magnitude, 
along the V axis. For any particular measurement point 
below G = Ge, £q is at most 16%. An average error for 
the whole fit, € = meang(€g), can also be inferred. To 
do so, we first have regularly under-sampled the experi- 
mental data onto 40 Jm wide bins, keeping there only 
the mean Gg and the mean log,,(V4). This way, and do- 
ing so separately for the two propagation branches (see 
appendix C, Fig. 15), no bias is introduced on = by the 
strong difference in measure density along the experi- 
mental (Va, Ga) curve (i.e., see Fig. 3). The thus derived 
overall fit error computes to = = 4%, below G = Ge. 
We discuss, in the next section, the fit beyond Ge and 


further discuss the accuracy of the inverted parameters 
in appendix B. 


B. On the fast crack velocity in PMMA 


Our simple sub-critical model hence matches most of 
the rupture dynamics of PMMA, from slow to fast ve- 
locities. In Fig. 3 however, an increase in velocity holds 
in the experimental data beyond G = G,, and is not 
properly accounted for. To highlight this mismatch, we 
display in Fig. 4 the fast branch with an optimised display 
scale. It has been shown [33] that, beyond a particular 
load, the global front velocity is impacted by the fracture 
instabilities that occur at high speed. Indeed, passed this 
threshold, fronts get more complex as micro-cracking oc- 
curs [32, 43, 44], that is, as micro-cracks form and propa- 
gate in the fracture plane ahead of the main front. Such 
micro-cracks are shown in Fig. 4. And, at an even higher 
load, micro-branching also comes into play, and aborted 
out-of-plane secondary cracks are observed [31, 58-60]. 
For the PMMA that is here studied, the micro-cracks 
were observed [32] at velocities above 165ms!, which 
approximately corresponds in the model to G > Ge. 
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FIG. 4. Zoom on the PMMA fast propagation branch pre- 
sented in Fig.3, and as per Eqs. (4) and (5). Beyond a load 
comparable to the modelled G. threshold, some micro-cracks 
start to nucleate, impacting the overall propagation velocity 
as explained by Guerra et al.[33]. The individual velocity 
of each micro-crack stays however constant at V = Vmicro- 
The validity of our single front model is limited passed this 
point, although it does predict a velocity plateau Vim, as per 
Eq. (10), and a velocity maximum Vmax, which are compara- 
ble to Vmicro. Inset (a): Fractography of the secondary micro- 
cracks on a postmortem fracture surface. White areas mark 
their nucleation centres. Inset (b): Atomic Force Microscopy 
of a nucleating cavity at the centre of a micro-crack. As pro- 
posed in section IV D, it could derive from the sublimation of 
localised bubbles around the main front, due to some intense 
thermal effects. 


Beyond this threshold, the apparent macroscopic speed 
of the front, V, increases with the micro-cracks grow- 
ing density, while the individual velocity of each micro- 
front, however, was inferred to stay constant [33], around 
Vinicro ~ 200ms~!, as illustrated in Fig. 4. 

Such a plateau in the propagation speed is somewhat 
consistent with our description (see Fig. 4). But this be- 
ing said, it is clear that our unique front model shows 
limitations as soon as fronts complexify. We can still 
push this discussion on the fast regime a bit further. 
A question of interest about the rupture of PMMA has 
been why the maximal observed crack velocity was sig- 
nificantly lower than the theoretical Rayleigh speed [32] 
(i.e., about 200ms™ rather than 880 mst). Equation 
(4) gives here some insight, as it does predict a plateau 
velocity Viim as the applied G gets very large. Indeed, 
besides preventing the crack advance at very low loads, 
the sub-critical healing processes significantly limit the 
fast growth rate, as the tip temperature is modelled to be 
high. More specifically, by inserting AT fast (7) in Eq. (4), 
and by looking at the high loads asymptotic regime of the 
healing term, we predict V to be limited by co = dC, 
the individual heat capacity of atom bounds: 


En p dé(1 + Gr/G) 
Vim ~ Vo|1 exp ( a 
(10) 
mco l 
~Vo}1 exp ( kp 5) 


In this expression, the crossed out terms are neglected 
in regard to the neighbouring ones. Note however that 
Eq. (10) is mainly illustrative, as the plateau it describes 
occurs in a domain where our single front model does not 
strictly apply. Note also that the value Vim ~ 100m st 
is smaller than the modelled maximum individual prop- 
agation velocity Vmax ~ 160ms!, which is obtained for 
G = Ge rather than for G — +00 (see Fig. 4). 


C. The detachment of Pressure Sensitive Adhesives 


We now pursue the comparison with the reported rup- 
ture of another material, acrylic based pressure sensitive 
adhesives (PSA), that typically happens when unrolling 
some office tape. In particular, the peeling dynamics of 
Scotch® 3M 600 rolls (composed of a polyolefin rigid 
backing coated with a layer of synthetic acrylic adhe- 
sive) has been thoroughly studied in the last decades 
(e.g., [19, 37, 61]); we here fit our model to two compatible 
(G, V) data sets that were published by Dalbe et al. [37] 
and by Barquins and Ciccotti [19]. These data sets are 
shown in Fig.5. Two stable modes of front detachment 
(i.e., a fast one and a slow one) are reported [19], similarly 
to those governing the rupture in PMMA. Additionally, 
some (unstable) stick-slip in the rupture dynamics is also 


observed [37] when peeling with an average velocity be- 
tween V ~ 15cms°! and V ~ 20ms!. 

Overlaying this experimental data, Fig.5 also displays 
a calibrated version of our model. The model parame- 
ters were inverted as follows, with a similar asymptotic 
analysis as what was done for PMMA. As no significant 
healing threshold displays at low velocity, we have only 
assumed that Ga is high enough to completely neglect 
the healing processes (i.e., the healing term in Eq. (4) is 
small if Ga is high). Of course, this absence of thresh- 
old, below which no forward propagation of the crack is 
observed, could also only indicate that (Ge — Gn)/2 < 0 
or that this value (i.e., illustrated on the PMMA data 
in Fig. 3) is less than the minimum energy release rate 
that was investigated in the tape experiments. We dis- 
cuss this particular point further in appendix D. We now 
invert the length £, which is, again, given by the slope of 
the slow phase and is here about 10nm. As no healing 
is now supposed to be at play, the nominal velocity Vo 
is given by the highest velocity records: Vo ~ 30ms"! as 
Vo is the maximum value then predicted by Eq. (4). Sat- 
isfyingly, this value compares well with the magnitude 
of a mechanical wave velocity in PSA, that is, 4/u/p, 
where u is, for instance, the shear modulus of the adhe- 
sive [62], 0.1 to 1 MPa, and p is its volumetric mass [63], 
about 10°kgm?. Next, from Eq. (8), the intercept of 
the slow branch with the ordinate (zero G) axis indicates 
Ge ~ 150Jm?. Rather logically, and with the inverted 
value of €, this again corresponds to a value of fracture 
energy Ue ~ 1eV. Note also that Ge is again higher than 
the transition load Ga at which a creeping front jumps 
to a fast regime. From this transition load, arising in the 
model from the temperature rise at low velocity (6), we 
also infer \/¢ to be in the order of 0.1 Jst mt K=. As 
the adhesive’s conductivity À lies in the same range [64], 
a consequent portion of G should be released into heat: 
o ~ 1. Of course, ¢ cannot be exactly one, as other dissi- 
pating processes than heat diffusion are likely to dissipate 
a part of G (see the discussion in section IV). According 
to our inversion however, this part ought to be small. Fi- 
nally, by varying / and by matching the coolest points of 
the fast phase, we estimate this parameter, which limits 
the highest tip temperature (i.e., Eq. (7)), to be in the 
nanometer range. This value, for the length scale of the 
heat production zone, is again rather consistent with the 
inverted magnitude of €, that is, the equivalent length 
scale for the mechanical shielding of the tip. Note also 
that both € and l are interestingly comparable to what 
was obtained for PMMA, and in the order of a polymeric 
entanglement density [57]. 

As shown in Fig.5 and with this set of parameters, the 
model accounts for most of the tape peeling dynamics. 
More quantitatively, for all the particular data points 
of the two stable phases, the fit error eq (as defined by 
Eq. (9)) is less than 20%. We also computed a mean fit 
error € = meang(€q) for the stable phases. To do so, and 
as done for PMMA, we first averaged the data points 
onto 10Jm wide bins, so that no densely populated 
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FIG. 5. Crack velocity V as a function of the energy release rate G as predicted by Eq. (4) and (5) (plain curve) and fitted to the 
tape experimental data [19, 37]. The unstable branch was not actually measured and the data points there are only averaged V 
versus G for a crack that undergoes stick-slip, in the given set-up, between the slow and the fast phase. The arrows indicate to 
which model parameters each part of the curve is mainly sensitive, and the main color scale specifies at which temperature the 
crack tip is modelled to be. The load Ga is an avalanche threshold beyond which peeling fronts can only propagate quickly and 
Ge is the modelled microscopic energy barrier for rupture. The inset shows the associated modelled temperature field around 
the front, at the onset of the fast to slow phase shift (G = 50J m?, V = 20ms"t). For readability, the color map is there 
different from the main one, and the circle corresponds to the tip of radius l, where the extra heat is emitted. AT of the main 


model curve is the value at the centre of the circle. 


part of the measured curve dominate the value of = (see 
appendix C, Fig. 16). We thus computed = = 5%. 

Note that, in comparison to the fast branch for the fail- 
ure of PMMA (i.e., as discussed in section IIIB), it would 
be of interest to know if the critical load G = Ge also ap- 
proximately corresponds to the apparition of some new 
rupture modes. Yet, the high velocity branch of the tape 
data is bound to relatively large uncertainties (the load- 
ing system of Barquins and Ciccotti [19] involved drop- 
ping weights from an elevated balcony, illustrating the 
challenges in fast peeling measurements), so that it does 
not allow a more thorough analysis. 


D. Parameter summary 


In Tab.I, we summarises all the parameters’s values, 
that we have inverted or supposed for the rupture of 
PMMA and PSA. The accuracy of these values is fur- 
ther discussed in appendix B. 


Parameter PMMA PSA Unit 
Vo 880 30 ms! 
Ge 1300 150 J m? 
Gh 650 z Jm? 
E 50 10 nm 
l 1 1 nm 
$ 0.2 ~1 
À 0.1 0.18 Jst mt kt 
C 1.5 1 MJ m Kt 
Ge/Ga 1.8 1.6 - 
U. 1 1 eV 
N 500 100 - 


TABLE I. Summary of all model parameters considered for 
the rupture of PMMA and PSA, as discussed in section III. 
A value do ~ 2A has been assumed in the derivation of 
these parameters. For completeness, the related quantities 
Ge/Ga = Uc/Ua, Ue = Gedÿ/(2Ë) and the shielding factor 
N = 2€/do are also specified. 


IV. DISCUSSIONS 


For two different polymeric materials, we thus have 
shown how a thermally activated fracture process, cou- 
pled with the dissipation and diffusion of heat, can sim- 
ply explain many features of the dynamics of both creep- 
ing and fast cracks, and the shifts from one state to the 
other. Such novel match, over seven to nine decades of 
propagation velocities and with only very simple physics 
considerations, could shade some new light on fracture 
mechanics, as thermal effects are often discarded. 


A. How hot is too hot for a crack tip? Some light 
from fractoluminescence 


To explain the fast propagation branch, we have no- 
tably predicted the front temperature to reach several 
thousands of degrees. Such high values are difficult to 
confirm experimentally, especially as they are to stand 
only on a few nanometers during short avalanches. There 
exist however, indirect hints toward the existence of an 
important temperature elevation in a variety of brittle 
materials fracturing at high speed. 

For instance, the analysis of some fracture roughness in 
cleaved quasi-crystals has revealed a damage zone of size 
anomalously large for this class of materials, and this 
was stated to result from a local temperature elevation 
of about 500K at the moving crack tip [65]. 

Several experimental works in glass and quartz [10-12] 
also managed to indirectly measure AT to indeed reach 
thousands of degrees, by characterising the photons emis- 
sion from the tips of some moving cracks and by com- 
paring it to the blackbody radiation theory [66]. In the 
case of tape, when peeling fast enough to be in the stick- 
slip regime, a blue tribo-radiation can similarly be ob- 
served [19, 67], and it was established that this radiation 
only occurs during the fast propagation phases of the cy- 
cle [19]. A direct example of such an emission is shown in 
Fig. 6, and its color could well correspond to the central 
wavelength Apeak associated, via Wien’s law [66], with a 
blackbody temperature compatible with our model: 


b 


= ~ 400 11 
TL AT nm, (11) 


Apeak = 
where b is Wien’s displacement constant ~ 0.0029mK 
and AT is about 7000 K at a load just passed the stick- 
slip threshold G = 90 Jm? (see Fig.5). The intensity of 
the observed light, which is visible in the dark but not 
under normal lightening, seems to also be consistent with 


the model. According to the Stefan—Boltzmann law [66], 
we indeed expect a radiated power in the order of 


P = s(To + AT)*hil ~ 1mW, (12) 


where s ~ 5.67x1078 W m? K4 is the Stefan-Boltzmann 
constant, and h the tape width (2cm) so that hl is 
the total area that significantly emits light. Note that 


such a power only accounts for a negligible part of 
the energy that is dissipated as the front advances, as 
P/(GVh) ~ 1074. For a human eye at a distance 
D ~ 10cm, it corresponds to a light luminance of about 
eP/(4r D?) ~ 1cdm?, using a blackbody luminous effi- 
cacy [68] e of 100 lumens per watt. With a pupil open- 
ing of about 10mm?, such a luminance is in the order 
of 10trolands (Td) [69], which does fit that of a flick- 
ering (i.e., the front has a stick-slip motion) radiation 
that is only visible in the dark, as those approximately 
range between 0.01 and 100 Td[69]. While the eye is 
persistent, a camera sensor of size S ~ 10mm?, placed 
at the same distance, would capture an averaged power 
yPS/(4r D?) ~ 100nW, where + ~ 0.1 is a typical ratio 
of time during which the front is in the fast phase com- 
pared to the total recording time, when peeling at a slow 
average velocity (i.e., ~ 15cm s`!) [37]. The magnitude of 
this power is interestingly close to the 10nW that were 
successfully measured by Camara et al. [67] for the lumi- 
nescence of another adhesive roll. 

For a given PMMA, Fuller et al. [70] also tried to quantify 
the temperature elevation around a quick fracture, both 
with the thermoluminescence technique and by using a 
liquid crystal coating on the matrix, whose color was 
thermosensitive [71]. For cracks propagating at 400 ms 
and faster, they measured heat efficiencies of about 2000 
Jm?, which is fairly compatible with the value we have 
derived for dG (a 400ms°! speed is obtained for G > 
4000 Jm? in Fig.3 and ¢ was inferred to be about 0.2). 
This experimental work [70] also estimated the instanta- 
neous temperature elevation of the fractures to be about 
500 K over a 0.5-um-thick area around the front. Such a 
thickness for the heat source was however acknowledged 
to be rather uncertain, as the measure sensibility for this 
parameter was limited. We remark that the same energy 
spread on the l ~ 10 nm thickness which we have here in- 
ferred would give a temperature rise of 104K and more, 
as predicted by our model (see Fig. 3). 

Truly, fractoluminescence could emanate from other 
mechanisms than some hot matter radiation. It was for 
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FIG. 6. Blue radiation emitted when quickly peeling tape 
beyond the stick-slip threshold (i.e., at an average velocity 
greater than 15 cms !, see Fig. 5). This picture was captured 
in the dark by a standard reflex camera (ISO: 25600, shutter 
speed: 1/2 s, focal length: 60 mm, aperture: f/4). The low 
shutter speed ensures that enough light enters the camera, but 
then covers many stick-slip cycles of the peeling dynamics [37]. 
Such fractoluminescence could be the mark of a very hot crack 
front [10-12] when unrolling tape. 


instance proposed [67, 72] that it partly arises from the 
molecules excitation of the fracture in situ air, by some 
electrical discharges between the two crack planes. Both 
these phenomena could surely coincide and, in any case, 
the light emission is an indication that some extreme and 
localised phenomena are at stake during fast failure. In 
that way, the thermodynamics model we propose holds 
some compatibility with that of Slepyan [2], where the 
abrupt advance of cracks derives from the emission of 
high frequency phonons, that excite atom bonds ahead 
of the tip, but that do not necessary thermalize. 

Some relatively recent atomistic simulations [73] seem 
nonetheless to confirm that the atoms at a moving front 
can undergo a significant heat. In a modelled graphene, 
Budarapu et al. [73] thus inferred a 200K temperature 
rise, over a 43 nm x 43nm area surrounding a running 
tip. This estimation is interestingly compatible with the 
thermal maps presented in Figs 3 and 5, for which the 
mean temperature is respectively 950 K and 350 K, when 
recomputed on a similar 1800 nm? surface upon the front. 
Note that atomistic simulations might naturally be more 
proper than our mesoscopic description, in particular be- 
cause the small scales (l) and high excitation frequen- 
cies (V/{) at play could call for more complicated mod- 
els [74, 75] than plain Fourier diffusion, Arrhenius growth 
or blackbody radiation. Yet, atomistic simulations are by 
nature far heavier to run, requiring an accurate descrip- 
tion of the atomic interactions onto femtosecond time 
steps. 


B. Is asimple model too simple? 


It is actually surprising that the proposed simple meso- 
scopic model can describe the propagation of cracks, 
when such a propagation, in reality, displays many com- 
plex phenomena. For instance, we have completely 
neglected the impact of crazing on the crack dynam- 
ics [36, 52], that is, the formation of defaults and fibrils 
at relatively large scales around the fracture front (i.e., a 
hundred of micrometers in PMMA and up to millimeters 
in PSA), while such large scale plasticity is often con- 
sidered to have a strong effect on the growth of cracks 
(e.g., [36, 52]). Yet, crazing is not incompatible with our 
thermal weakening model, which only states that a signif- 
icant part of the mechanical energy should be dissipated 
far closer to the crack front (i.e., over a few nanome- 
ters), and that this very local dissipation should be that 
of a first effect on the crack dynamics. In this descrip- 
tion, crazing is then a consequence of the front progres- 
sion rather than its main cause. In a similar way, many 
other known failure phenomena, such as the emission of 
mechanical waves during rupture [50], complicated creep 
laws from the corrosive interactions between the frac- 
ture fluid and the fracture tip (e.g., [46], chpt. 5.4), or 
the complexification of fronts at high propagation veloc- 
ities [44], are not directly encompassed by Eqs. (4) and 
(5), but are not in conflict with the model either. 
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The simplicity of the model can actually be considered 
as one of its strength, as the physics that it describes 
could apply to many different materials and not only to 
polymers. Accurately testing this idea would however 
require the full (G, V) curves of more materials, and 
those are often not trivial to obtain experimentally at 
all velocities. Such experimental work could yet be re- 
warding, as we have here shown that matching the model 
to some (G, V) curves can give some valuable insights 
on the rupture of matter. Our quantification for each 
model parameter stays however rather approximate, and 
we have mainly derived their orders of magnitude. We 
have, in particular, assumed that they were all constant 
for a given material, while most could be velocity or tem- 
perature dependent [18, 24, 76]. For instance, the fact 
that PSA exhibits a larger scale viscous behavior (i.e., 
including fibrillation and heating over millimeters around 
the tip) at lower velocity [36] could indicate that the heat 
production size | decreases with the crack speed in this 
medium. It is especially known that the elastic moduli 
in PSA are strongly temperature dependent [62], and this 
was actually proposed by Maugis [21] and Carbone and 
Persson [24] as the driving cause for failure instability in 
rubber-like materials. We have, besides, considered both 
PMMA and PSA as homogeneously tough while Ge is 
bound to present some quenched disorder. While such 
heterogeneities should not affect the stable propagation 
branches, as long as G and V are then understood quanti- 
ties which are averaged over a few Ge correlation lengths, 
it could be of importance for the accuracy of the loads 
at which the phase transitions occur [28], as slow cracks 
shall preferentially avalanche on weaker zones and fast 
cracks stop on stronger locations. In the case of PSA, we 
have furthermore considered that peeling was a cohesive 
process (i.e., that it occurs inside the adhesive), while a 
bi-materials interfacial model would be more appropri- 
ate, as the crack essentially propagates at the interface 
between the substrate and the glue [77]. 

These numerous limitations being stated, the parameters 
we have inverted are nonetheless in rather satisfying or- 
ders of magnitude, confirming the physical relevance of 
the model. Indeed, the intrinsic fracture energy in both 
materials Ue = dÿG./(2£) is comparable to one electron- 
volt, which is typical for an energy that bonds atoms 
(e.g., see appx. E in Ref.[6]). Because our proposed 
description is statistical, one should remember that Ue 
is a mean material feature, for a rupture process that 
is made of several types of bond breaking. As a rough 
example, Ua ~ 1eV may indicate that the crack consum- 
mates in average three weak links (such as hydrogen or 
Van des Waals bonds of respective energies [6] ~ 0.1 and 
~ 0.01 eV) for every stronger connection that snaps (say, 
one C-C link of an acrylic chain, of covalence energy [6] 
~ 4eV). The nanometric scale / for the heat generation 
may well correspond to the typical entanglement density 
in polymers [57] (the density of polymeric chains cross- 
ing points in the matrix), below which atoms have more 
freedom to vibrate, and which is known to affect some 


rupture properties (e.g., [52, 57]). It is also coherent that 
the generation of heat was inferred to occur over a length 
scale comparable to €, the equivalent radius describing 
the energy shielding of the tip. We have indeed derived 
that the former is a strong cause for the latter, as the 
heat efficiency ¢ was inverted to be non negligible (i.e., 
@~0.2—1). 


C. Tip stress and front shielding 


A nanometric scale (i.e., comparable to € or l) has been 
noteworthily observed in the rupture of other materials. 
One example is the length scale of a light radiating (and 
hence likely thermal) zone around running fracture tips 
in glass [42]. In carbonate rocks, it is also the typical size 
of some observed nanograins that form along sliding seis- 
mic fault planes [78]. Such a nano-damage explains the 
glossy and reflective aspect displayed by some faults (of- 
ten referred to as fault mirrors), as their typical surface 
roughness is then comparable to the wavelengths of visi- 
ble light. The origin of this damage, however, is debated 
as, below 1 um, plasticity is expected to dominate over 
brittleness in this material and asperities should hence 
deform rather than break. Noteworthily, some intense 
thermal effects, arising from the frictional heat, such as 
some fast melting and cooling or the thermal decomposi- 
tion of carbonates, were proposed to solve this apparent 
paradox [79]. 

Similarly, for the materials that we have here studied, 
the usual predictions for the size of the shielding pro- 
cess zones are far larger than €. In PMMA, for instance, 
it is in the order of macro ~ GE/ø} ~ 200m, where 
oy ~ 100 MPa is the tensile yield stress of the bulk poly- 
mer and E ~ 3GPa its Young modulus [55]. However, 
in that description, oy is a stress that is averaged over 
a macroscopic sample, and is likely not representative 
of the actual energy density around the defaults of this 
sample. It was notably reported that a Dugdale [80] like 
cohesion model (i.e., ø is homogeneously equal to o, in 
a process zone of radius macro), poorly accounts for fast 
rupture in PMMA [81]. Naturally, macro is still to bear 
some significance, in particular as a characteristic length 
scale for crazing in acrylic glass [52], where a portion 
within (1 — ¢) of the release rate G is to be dissipated, 
either by the creation of dislocations [51], the emission of 
waves [42, 50] or residual thermal effects. But € was in- 
verted as an equivalent size, only defined by 2£/do = N 
with N the damping of the tip potential energy U due 
to the energy dissipation. We solely inverted N to be 
around 100 and 500 for respectively PSA and PMMA, 
and many links might well snap and heal far away from 
the tip, allowing for crazing. 

Still, most of the rupture is likely to occur very close to 
the front where the stress is to be the highest. We can 
estimate such a stress at the tip by considering a simpli- 
fied expression for the elastic energy stored in rupturing 
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(13) 
which, with Eq. (3), is equivalent to the well known form 
for the limitation of an otherwise divergent stress at the 
tip of cracks, predicted by the general elasticity theory 


(e.g., [46]): 
o~ 4] —. (14) 


In the case of PMMA, such a computed stress is as high 
as 7 GPa, and we thus predict a high atomic strain o/E of 
about 200% at the onset to fast rupture (i.e., for U equal 
to U./1.8 as per Tab. I). Such a strain shall be likely at 
a fracture tip for the strong intermolecular deformation 
immediately before failure. Of course, the simply linear 
elastic Eq. (13) is unlikely to be valid at 200% strain, 
and we also considered describing U with a Morse poten- 
tial [82], that is 


zafi ep soa w) (15) 


which, at the onset of fast rupture, predicts a strain 
(d — do)/do ~ 400%. While this dual-particles poten- 
tial stays, by nature, a strong approximation in the com- 
plex rupture of a polymer, it is worth reminding that the 
model which we have introduced does not rely on a par- 
ticular shape of the inter-atomic potentials (i.e., neither 
on Eq. (13) or on Eq.(15)), but only on their average 
dissociation energy Ue. 


D. Front complexification 


Overall, our derivation of é < macro only suggests 
that process zones are heterogeneous objects, dissipat- 
ing a higher density of energy in their centre than at 
their periphery. In particular, it was shown that a few 
tens of micrometers (i.e., a portion of macro) is a typi- 
cal distance at which the secondary micro-cracks nucle- 
ate from the main front in PMMA [33] and, as shown in 
Fig. 4, the imaging of some postmortem rupture surfaces 
reveals that these micro-cracks initially grow from iso- 
lated spherical cavities at their centre, of radius about 
300nm. We here propose that such cavities could cor- 
respond to bubbles, forming by sublimation [83] on weak 
locations of the process zone, and leading to some micro- 
fractures once having grown to a critical size. While 
remaining to be confirmed, such a sublimation process 
would definitely require some local but very high temper- 
atures in the crazing area. Indeed, to nucleate ahead of 
the main front, the observed cavities have to form during 
less than Emacro/V ~ 1 us, and the pyrolysis of PMMA to 
methyl methacrylate (MMA) only reaches such a reaction 


rate at temperatures T, that are beyond 1000 kelvins [83]. 
In return, and assuming that the ideal gas law approx- 
imately applies (e.g., [17], chpt. 4), some bubbles form- 
ing at this temperature would hold an internal pressure 
pRT,/M, where M = 0.1 kg molt is the MMA molecular 
mass [84], p = 1200kgm is the volumetric mass of the 
solid PMMA [55] and R is the ideal gas constant. This 
value computes to at least 100 MPa, which is comparable 
to the surrounding bulk compressive strength [55]. The 
evolution from pressurised pores to propagating micro- 
cracks would then be coherent. 

Thus, in addition to explaining, as shown in this work, 
the first order dynamics of singular fronts, concentrated 
thermal processes could also be responsible for their com- 
plexification at high propagation velocities. In the case of 
acrylic glass, we have notably inferred (see section III B) 
that the appearance of the secondary fronts approxi- 
mately coincides with energy release rates that are close 
to the (modelled) intrinsic barrier Ge. This concomitance 
could be explained by the need for new dissipation pro- 
cesses, when cracks propagate over-critically (G > Ge) 
so that some extra energy is brought to the rupture sys- 
tem. Such an idea is notably re-enforced by the fact that 
the density of nucleated micro-cracks was inferred to be 
proportional to a value comparable to G — Ge, as shown 
by Guerra et al. [33]. 


V. CONCLUSION AND PERSPECTIVES 


We presented a new and general model for the kinetics 
of cracks. The main physical elements that were intro- 
duced in this model are, only, a sub-critical (Arrhenius- 
like) growth rate and the dissipation and diffusion of heat 
around fracture tips (where the applied mechanical stress 
is concentrated), an immediate consequence of the latter 
being the possibility for crack fronts to reach thousands 
of degrees temperatures. Interestingly, these different ele- 
ments have, separately, long been considered or observed 
in the physics of rupture (e.g.,[7, 18, 39, 85]), but had 
not previously been combined for comparison with some 
experimental data. In doing so, we here showed that the 
rupture of two materials (namely, PMMA and PSA) can 
be quantitatively reproduced over many decades of prop- 
agation velocities, from slow creep regime to fast propa- 
gation. 

Thus, we inferred that the propagation of a crack can 
be sub-critical, even at velocities approaching that of the 
mechanical waves in the surrounding matrix, due to its 
potentially very high tip temperature. We also suggested 
that the microscopic healing process around a fracture 
front can significantly constrain the fast velocity regime, 
from the strong thermal activation at such temperatures, 
while it is often considered that healing is only relevant 
for very slow cracks. The existence of thousands of de- 
gree temperatures is actually supported by many exper- 
imental works that study the visible fractoluminescence 
of fast fronts [40—42, 70]. In some instance [73], it has also 
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been modelled by some atomistic simulations, and we ad- 
ditionally showed, in the present work, the existence of 
bubble forming in the process zones of cracks in PMMA. 
We proposed that these bubbles could well originate from 
some local sublimation of the polymer near crack tips. As 
they are located at the nucleation centres of secondary 
fracture fronts, we also suggested that the complexifica- 
tion of cracks at high velocities could derive, as the rest 
of the propagation dynamics, from some thermally acti- 
vated processes. Finally, for the two materials that we 
have studied, we have inferred that the mechanical stress 
around cracks remains an increasing quantity inside the 
process zones up to a few nanometers from the tip. Such 
a nanometric scale matches the typical size over which 
the heat was inferred to be generated, making thermal 
dissipation the likely main process that shields rupture 
fronts from mechanical failure. 

Noteworthily, the proposed model, and its ability to 
explain some actual crack dynamics, stresses the im- 
portance of the heat conductivity of materials on their 
macroscopic strength. A high conductivity indeed al- 
lows to evacuate the extra internal energy away from 
the fronts, thus delaying any thermal weakening. As a 
general statement, many strong materials happen to be 
good conductors, such as metals, graphene [86] or spider 
silk [87]. For the latter, it was in particular shown that, 
contrarily to most materials, its conductivity actually in- 
creases with deformation [87], which could well be a natu- 
ral defence mechanism for the stability of arachnid webs. 
Designing human-made solid matrices that can replicate 
such a behavior on the molecular scale could then be- 
come a new important target of material sciences. 
Finally, we suggest that most of the physics that we have 
introduced to study mode I fractures shall also be valid 
for mixed-mode fracturing as well as for solid friction. 
The latter is actually suspected to hold some non negli- 
gible, thermal related, weakening mechanisms (e.g., [88]), 
which could notably be a key in geophysics in under- 
standing the stability of seismic faults. Such mechanisms 
might be diverse, and may include the thermal pressuri- 
sation of fault fluids [89, 90] or some changes in the fault 
planes minerals phase (i.e., such as melting or thermal 
decomposition) [91]. We propose that they could also be 
related to a thermally boosted sub-critical slip, in the 
sense of statistical physics and similarly to the model we 
have here developed. 
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Appendix A: Method for the measure of crack 
velocity versus energy release rate in PMMA 


As part of the PMMA data was not published before 

(i.e., the slow propagation branch), we describe, in this 
section, the method that was used to acquire it. 
Wedge splitting fracture tests are used to measure both 
the slow and fast V(G) branches in PMMA [32, 54], 
whose geometry is shown in Fig. 7. Rectangular plates of 
size 140mm x 125mm x 15mm are first machined from a 
plate of moulded PMMA (Perspex®). A 25mm x 25mm 
notch is subsequently cut out on one of the two lateral 
edges and a 8-mm-long 800-um-thick groove is finally in- 
troduced in the middle of the notch with a diamond saw. 
To grow slow cracks, an additional seed crack (~ 2 mm- 
long) is added at the end of the groove via a razor 
blade. This crack is loaded in tension by pushing a steel 
wedge (semi-angle of 15°) in the notch. Two steel blocks 
equipped with rollers are placed in between the wedge 
and the specimen notch to limit the parasitic mechani- 
cal dissipation through plastic deformations or friction at 
loading contacts. As a result, the vicinity of the crack tip 
can be assumed to be the sole dissipation source for me- 
chanical energy in the system. The wedge speed is first 
set to 1.6yms!. The force F, applied by the wedge to 
the specimen, increases linearly with time up to a point 
F above which the seed crack starts to propagate. Above 
this point, F decreases with time. We let the crack propa- 
gate over a distance of about 10mm. This ensures repro- 
ducible initial conditions with a long-enough well-defined 
sharp seed crack. The specimen is then unloaded (un- 
loading wedge speed: 16yms). The specimen is then 
loaded again at a constant prescribed wedge speed Visage: 
which has been varied from 1.6 ms"! to 1.2mms!. 
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FIG. 7. Schematic of the experimental set-up used to measure 
the crack energy release rate G and its corresponding propa- 
gation velocity V(t) in PMMA. See Refs. [32, 54] for details. 
The hole is used to store some potential energy in the PMMA 
sheet for fast propagation experiments and is replaced by only 
a seed crack for slow ones. The dashed horizontal lines rep- 
resent conductive metallic lines deposited onto the sample to 
measure the fast crack velocity with an oscilloscope. 


During each fracture test, the force F(t) is monitored in 
real-time via a cell force mounted on the system (S-type 
Vishay load cell). A camera (USB2 uEye from IDS) is 
also used to image crack propagation at the specimen 
surface (space and time accuracy of 125 um and 0.1s). 
A coarse approximation of the crack speed can be ob- 
tained by differentiating the position of the crack tip ob- 
served on the successive images. However, a more ac- 
curate signal V(t) is obtained from the force signal (see 
Ref. [92] for details on the method). Indeed, in a lin- 
ear elastic isotropic material like PMMA, the specimen 
stiffness k(t) = F(t)/(Vweaget) is a continuous decreasing 
function of the crack length, c(t), that is set by the spec- 
imen geometry only. This function has been obtained 
using finite element calculations on the exact experimen- 
tal geometry (Cast3M software, 2D simulation assuming 
plane stress conditions); it was checked that the obtained 
k versus c curve coincides with the experimental curves 
obtained by plotting k(t) as a function of the crack length 
measured by the camera. The idea is then to use this 
curve k(c), and the corresponding inverse function k~?, 
to infer the time evolution of crack length c(t) from the 
signal F(t): c(t) = k71[F(t)/(Vweaget)|. Time deriva- 
tion of the so-obtained c(t) provides a signal V(t) about 
50 times less noisy than that directly obtained from the 
camera images. The knowledge of c(t) and F(t) also al- 
lows determining the time evolution of the energy release 
rate, G(t). Indeed, the total amount of mechanical en- 
ergy provided to the specimen is F?(t)/[2k(c(t))]. Dif- 


ferentiating this stored energy with respect to c directly 
provides G(t). The slow branch of Fig.3 then provides 
the observed V(t) as a function of G(t). The results from 
twelve fracture experiments are gathered in this branch 
and differ by their Vweage value. 

To grow fast cracks and measure V(G) in the fast sta- 
ble phase, the seed crack has been replaced by a hole 
of tunable radius (1 to 4mm) drilled at the end of the 
groove [32]. This delays fracture and increases the po- 
tential energy stored in the specimen at the initiation of 
crack growth. The time evolution of V(t) is measured by 
monitoring, via an oscilloscope, the successive rupture of 
parallel 500-um-large metallic lines (chromium/gold) de- 
posited on the surface. That of the stress intensity factor 
K is obtained via finite element analysis (see Ref. [32] for 
details). The time evolution of the mechanical energy re- 
lease rate is then deduced: G = K?/E where the Young 
modulus ŒE in the studied PMMA have been measured 
to be E = 2.8GPa. The fast branch of Fig.3 then pro- 
vides the observed V(t) as a function of G(t). The re- 
sults from five fracture experiments are gathered in this 
branch; they differ by the amount of stored elastic energy 
at crack growth initiation. 


Appendix B: Parameters sensitivity 


We here show, on the PMMA data, how varying the 
model parameters around their inferred values impacts 
the model fit, thus giving the reader a better feeling for 
their individual effect and sensitivity. In each of the fig- 
ures 8 to 14, a unique parameter of the model varies 
while the others are kept to the exact values used for 
the fit presented in Fig.3: € = 56nm, Vo = 880ms'!, 
Ge = 1275Jm?, Gah = 650Jm?, à = 20%, À = 
0.18Jst mt Kt, C = 1.5 x 106 Jm? K1, l = 1nm and 
To = 296 K. These seven plots show the fits up to the ap- 
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FIG. 8. Effect of varying the nominal velocity, Vo, on the 
fit to the PMMA data. The propagation velocity is roughly 
proportional to Vo, but also modifies the positions of the phase 
transitions. 
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FIG. 9. Effect of varying the breaking energy barrier, Ge, on 
the fit to the PMMA data. At a given load G, the higher 
Ge, the slower the crack. The transitions between the three 
propagation modes (fast, slow, and dominated by healing) 
are also affected: a medium with a stronger barrier needs a 
heavier load to transit to a weaker state. 
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FIG. 10. Effect of varying the healing energy barrier, Gn, on 
the fit to the PMMA data. A crack that heals more easily 
needs a higher load to actually propagate forward or to stay 
in the high velocity regime. 


parition of the secondary micro-cracks (see section IIIB 
and Fig. 4), after which the model does not apply as such. 
Naturally, some care should be taken when interpreting 
the inverted parameters (i.e., €, Ge, Gh, @ and l) beyond 
their actual orders of magnitude. For instance, € and Ge 
were fitted by a linear regression (i.e., Eq. (8)) on the data 
which lies between G = 350 and G = 700 J m? in Fig. 3, 
and the above values (£ = 56nm and Ge = 1275 Jm?) 
were obtained with a coefficient of determination R? 
equal to 0.85. Allowing R? to drop down to 0.75 during 
this fit gives € in a 30 to 80 nm range and Ge between 950 
and 1500 Jm?. G'n being directly deduced from Ge (and 
from the vertical asymptote at low velocity in Fig. 3), 
the uncertainty on its value is comparable to that of Ge, 
that is, a few hundreds of joules per square meter. Let us 
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FIG. 11. Effect of varying the stress cut-off scale €, on the 
fit to the PMMA data. € mainly controls the slope and the 
intercept of the low velocity branch. A small change in € 
significantly modifies this branch as well as the threshold to 
the fast regime. 
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FIG. 12. Effect of varying the ratio of energy converted to 
heat, ¢, on the fit to the PMMA data. The maximum ve- 
locity increases with @ as the tip temperature is higher. The 
threshold from the slow to the fast branch (i.e., the (Ga, Va) 
point) shifts towards a lower G as a lighter load is required 
for the temperature to significantly deviate from To. 


now assess the accuracy of ¢. In the model, this parame- 
ters mainly controls which is the fastest point of the slow 
velocity branch (e.g., as shown in Fig. 12), after which 
cracks have to avalanche [28]. We then compare the ex- 
perimental value for this particular point (obtained at 
G = Ga and V = Vj, see Fig 12) to the model prediction 
of the same point. We quantify the error there as the 
euclidean distance between these points, in the sense of 
Eq. (9). Such a relative, unitless, error minimizes to 10% 
for @ = 0.25 and, should we allow it to rise up to 30%, 
we obtain ¢ to be between 0.15 and 0.30. Finally, let us 
assess the accuracy of the inversion for the length scale of 
the heat production zone l. We vary | and compute the 
same relative euclidean error in average over the fast ve- 
locity branch (i.e., the location where the model is mostly 
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FIG. 13. Effect of varying the thermal conductivity, A, on the 
fit to the PMMA data. With a higher À, the heat is better 
evacuated: the slow to fast branch threshold shifts towards 
higher G and V. The fast regime is not very sensitive to A, 
as AT is there constrained by I. 
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FIG. 14. Effect of varying the heat production zone radius, 

l, on the fit to the PMMA data. l mainly impacts the plot 

curvature on the high velocity branch. No effect is observed on 


the slow branch, as the thermal elevation there is constrained 
by the diffusion skin depth rather than by l (see Eq. (6)). 


sensitive to l, see Fig. 14). This error now minimizes to 
5% for l = 1 nm and, letting it reach 30%, we infer / to 
lie in a 0.1 nm to 2 nm range. 

These uncertainties in the parameter inversion are some- 
what high, but we here quantify an atomic scale process 
based on macroscopic measurements, so that this is not 
particularly surprising. One also needs to add up the ex- 
perimental inaccuracy for V and G (see the data spread 
in Fig. 3), as well as the limitations of our very first order 
physical model, as discussed in section IV B. Still, overall, 
the data is well explained over eight decades of velocities 
and with parameters that are in physically reasonable 
orders of magnitude. 


Appendix C: Data binning to compute a mean fit 
error 


To compute a mean fit error € = meang(€a), where £a 
is defined by Eq. (9), we first binned the PMMA and PSA 
data points onto coarse bins using a running average on 
both stable branches, as explained in the core text (see 
section III). This was done to avoid the densely populated 
parts of the data sets to dominate on the value of the 
inferred €. Figures 15 and 16 show the result of this data 
decimation. 
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FIG. 15. Under sampled PMMA data (blue circles), using a 
running average on 40 Jm? bins. The dots are the original 
data points and the black line the fitted model. With this 
decimated data set, the mean fitting error is = = 4%. The 
data is only shown to the onset of micro-cracking (i.e., see 
section III B), beyond which the model does not apply as such. 
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FIG. 16. Under sampled PSA data (blue circles), using a 
running average on 10 Jm? bins. The dots are the original 
data points and the black line the fitted model. With this 
decimated data set, the mean fitting error is € = 5%. 
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Appendix D: Healing processes in tape 


We considered the healing processes to be negligible 
in order to describe the dynamics of unrolling tape, 
as no low velocity constant G asymptote arising from 
crack healing displays in the (G, V) data (i.e., in Fig. 5). 
Such an absence would, however, also happen if Ge was 
to be smaller than Gp, as the asymptote is obtained 
for (Ge — Ch)/2. Thus, the healing energy barrier 
could still be comparable to the breaking one, and so 
still significantly impact the high velocity propagation 
branch, when the crack tip is hot enough for healing to 
be non negligible (as predicted by Eq. (4)). Of course, 
an accurate quantification of this effect suffers from 
the absence of the asymptote as it is the only good 
constraint for Gh. Figure 17 shows for instance a model 
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FIG. 17. Fit of the Scotch® 3M 600 data [19, 37] with a model 
including healing processes. The unstable (middle) branch of 
the model should not necessarily match the data point which 
are averaged G and V values for a front that stick-slips. 
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FIG. 18. Fit of the peeling data for another roller tape: 


Scotch® 3M 602 data [93]. The lack of linearity at low ve- 
locity calls for healing processes in our model. Note also the 
curvature on the lower end of the high velocity branch, not 
present in the other data sets that we show but rather com- 
patible with our proposed model. 


not disregarding healing, and compares it with the tape 


data. The match is there improved compared to the fit 
presented in Fig. 5 as we have now an additional degree of 
freedom. The fit parameters in this figure are as follow: 
€=9nm, Vo = 70ms!, Ge = 154J m?, Gr = 200 J m?, 
@é~1,\=02Js'm'!K!, C = 10° Jw? Kt, l =1nm 
and 76 = 296 K. 

Note that Barquins et al. [93], who released part of 
the data presented in Fig.17, also provided similar 
measurements for another type of roller tape, Scotch® 
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